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Abstract 

From a covariant Hamiltonian formulation, using symplectic ideas, we obtain covari- 
ant quasilocal energy-momentum boundary expressions for general gravity theories. 
£vq \ The expressions depend upon which variables are fixed on the boundary, a reference 

configuration and a displacement vector field. We consider applications to Einstein's 
theory, black hole thermodynamics and alternate spinor expressions. 
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1 Introduction 



The source of gravity is the energy-momentum density of all other physical 
fields. For the gravitational field itself, however, energy-momentum is not so 
simply described. Although geometries with suitable asymptotic regions have 
a well defined total energy-momentum, the equivalence principle precludes any 
proper local density. Hence, the idea of quasilocal quantities for gravity has 
been advocated and there have been many recent proposals. Here we present 
a set of quasilocal energy-momentum expressions for quite general geometric 
gravity theories. Distinguishing features of our Hamiltonian based expressions 
are that they are covariant and that the different possible expressions are 
associated with the choice of variables to be held fixed on the boundary and 
a symplectic boundary variational structure. 
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2 The covariant Hamiltonian 



We consider general theories of dynamic geometry. The possible geometric 
potentials are the coframe d a and connection u a p one-forms and the met- 
ric coefficients g^. The corresponding field strengths are the torsion Q := 
Dd a = aid" + u a j3 A d 13 and curvature Vt a p := duu a p + u a 1 A uo 1 ' p 2-forms and 
the nonmetricity one-form Q^ v := Dg^ v = dg^ u —ou a ^g ai/ — 0J a v g^- In place of 
the usual "second order" Lagrangian 4-form C = C(g, ■#, u, Q, 6, Q), we intro- 
duce covariant canonically conjugate momenta forms to obtain a "first order" 
Lagrangian 4-form 

C := Q^Ait^ + Q a Ar a + Sl a ApJ - A(g, ■&; tt, r, p), (1) 



which yields "first order" field equations via independent variations of the 
fields g, $, uj and momenta n, r, p. 

For any fixed slicing of space-time by t=constant surfaces X along with a 
connecting vector field N, the decomposition of the first order Lagrangian 4- 
form according to the general pattern £ = dt A i^C- = dt A (£n^P A p — Ti(N)) 
identifies the covariant Hamiltonian 3-form (or density) [1], i.e., the generator 
of evolution along N: 



H{N)=i N A + g^Ai N n^ - i N $ a DT a - Q a A i^r a - n a p Ai N pJ 

- t N u%(D P J - g a y» - g m ^ + ¥ A r a ) + dB(N), (2) 

where B(N) := iN$ a T a + i^u a ppd 3 '. 



3 Covariant quasilocal expressions 



Noether's theorem for a translation applied to C yields the differential iden- 
tity d7i(N) = (terms oc field equations) and the algebraic identity 7i(N) = 
(terms oc field equations) +dB(N). Consequently, the Hamiltonian j s H(N), 
the integral of the density over a finite spatial region, has a conserved value on 
a solution given via Stokes theorem by J> aE B(N). This boundary integral gives 
the quasilocal quantities (energy etc.); the limiting value at infinity should be 
a total conserved quantity for N asymptotically Killing. 

However, the value of the Hamiltonian is not yet firmly fixed. As with other 
Noether currents we can add any exact differential thereby modifying the 
expression given above for the Hamiltonian boundary term. Moreover, such 



an adjustment only affects the boundary term in the variation of the Hamil- 
tonian and thus does not affect the basic physics: the dynamic equations. 
The Hamiltonian formulation, however, includes a further principle which 
fixes the possible forms of B. The proper form is identified by consider- 
ing the variation of the Hamiltonian density: 57i(N) = (field eq. terms) 
+di N (5g llu ii f " u + 5§ a A r a + 8u a p A pj). The total differential term in STi(N) 
produces a boundary integral which reflects both the choice of variables to be 
held fixed ("control" variables) and the symplectic structure [2]. The boundary 
term in 57i(N) vanishes (as it should) for fixed control variables on a finite 
boundary. But for a boundary at infinity the limit of the asymptotic fall offs 
generally gives a nonvanishing results (e.g., for Einstein's theory) then the 
above B(N) expression needs adjustment [3]. 

An acceptable boundary expression for general metric compatible gravity the- 
ories was found [1] and soon improved by Hecht [4]. From his work we saw 
more possibilities. The expressions for the Hamiltonian boundary term which 
are 4-covariant and yield a 4-covariant symplectic structure for the boundary 
term in S7i(N) were identified [5]. There are two types: 

By(N) := i N tp A Ap - qAtp A i N p, (3) 

B p (N) :=i N $AAp- qAtp A i N p, (4) 

depending upon whether the configuration field k-form tp or its conjugate mo- 
menta p is "controlled" on the boundary [2] (i.e., either Dirichlet or Neumann 
boundary conditions). Here tp and p are the values in a reference configuration, 
Atp — tp — tp, Ap = p — p and <; = (— l) k . Thus for the geometric fields 



, -Ag^ijrK^ I I i N $ a Ar a + Ad a A i N r a 
-Ag^i N 7t^ J \ i N hAr a + A$ a A i N r a 

i N uj a pAp/ + Auo a p A i N pJ 1 
i N u a (3A P J + Au a p A i NP J J 

where for each bracket the upper (lower) line is to be used if the field (mo- 
mentum) is controlled. Hence, as in thermodynamics, there are several kinds 
of "energy" , each corresponds to the work done in a different (ideal) physical 
process [2]. Our quasilocal boundary expressions are covariant — aside from 
the manifestly non-covariant explicit connection terms in B. These terms in- 
clude a real physical effect plus an unphysical dynamical reference frame effect. 
These effects can be separated using the identity 

(i N u a p)$ p = i N Q a + DN a - £ N $ a . (6) 



Via this identity the boundary term contains a time derivative of certain frame 
components. Such terms have been noted previously in Einstein's theory [6]. 

The variation of the Hamiltonian, in addition to the field equation terms, now 
includes for each variable one of the boundary terms 

cHn(5(p A Ap), or di N ( — Aip A 8p), (7) 



which reflect the symplectic structure and the control mode. Specifically, for 
the geometric variables, the total differential term in S7i(N) is of the form 
(UnC where 

' Sg^Air^ 1 | 5ti a A Ar a 1 J Scv a p A ApJ 
-Ag^SiT^ | [ -A$ a A5r a \ \ -Au a p A SpJ 



here again the upper (lower) line in each bracket corresponds to controlling 
the field (momentum). Our quasilocal expressions are uniquely determined by 
the Hamiltonian variation (8) and the requirement that all of the quasilocal 
quantities vanish when the fields have the reference configuration values. 

This general formalism readily specializes to coordinate or orthonormal frames 
and Riemannian, Riemann-Cartan or teleparallel geometry. For General Rela- 
tivity and the Poincare Gauge Theory (asymptotically flat or constant curva- 
ture) our expressions reduce to known ones [7,8] and give the total quantities 
at spatial and future null infinity [9]. Non- vanishing reference configurations 
(e.g., Minkowski or de Sitter metric, frames, connections) play an essential role 
in obtaining these total values. Our quasilocal expressions likewise depend on 
a reference configuration. A reasonable choice for the reference configuration 
is to embed the spatial surface and its boundary into a Minkowski space [10]; 
more generally one could use (anti) de Sitter space, a homogeneous cosmol- 
ogy, a Schwarzschild solution, etc. The evolution vector field can be selected 
to correspond to a Killing field of the reference configuration to obtain the 
quasilocal energy- momentum (and angular momentum). As in Ref. [10], alter- 
nate choices of the evolution vector field can be used to distinguish between 
quasilocal quantities and conserved charges, mass and energy, etc. 



4 Einstein's theory and a spherically symmetric example 



For Einstein's theory our quasilocal expressions reduce to only two possibilities 
depending upon the choice of boundary conditions — the frame (i.e., the metric 
or intrinsic geometry) or the connection (the extrinsic geometry) could be held 



fixed on the boundary. These Dirichlet and Neumann quasilocal expressions 
are 

B<> = i N u a ^AeJ + Au> a A i N eJ, (9) 

B u = i N uj a AeJ + Acu a p A i N °eJ, (10) 

with the corresponding boundary terms in the Hamiltonian variation 

87U(N)mdi N (-Aw a i3A6e a (} ), (11) 

SHUN) « di N (Su a p A Ae/), (12) 

respectively, where e a/3 = *($ Q A ^ /3 ). The expression (9) differs slightly from 
one due to Katz [11]. It is also similar to the famous expression of Brown and 
York [10]. Significant differences from the latter are that (i) our expression is 
covariant, (ii) our expression is not specialized to the timelike boundary being 
orthogonal to the spacelike hypersurface, (iii) our expression includes interac- 
tion terms between the physical system and the reference configuration, and 
(iv) we have relaxed the seemingly natural choice lapse = 1 as it prevents the 
quasilocal energy from attaining to the total energy in the limit for asymp- 
totically anti-de Sitter solutions. (There is now a new work [12] which also 
considers quasilocal quantities in asymptotically non-flat spaces.) 

Consider the static spherically symmetric metric ds 2 = — e 2 *o?£ 2 + e" 2 ^dr 2 + 
r 2 (d$ 2 + sin 2 6d(j) 2 ) where e 2 * = 1 — 2m jr + \r 2 and the reference con- 
figuration is given by m = 0, i.e., Minkowski or (anti) de Sitter. For the 
quasilocal energy within a centered sphere we find E$ = ra(e*° — e*) and 
E u = rae*~*°(e*°— e*). An appropriate choice of lapse is a = e* = VI + Xr 2 
which corresponds to the reference configuration timelike Killing vector. For 
this choice both expressions give m for the total energy, while the energy 
within the horizon is 2m for E$ and oo for E w . The quasilocal energy exterior 
to the horizon is negative in both cases. Note that there is no simple rela- 
tionship between the quasilocal energies of Schwarzschild anti-de Sitter space 
referenced to Minkowski space, Schwarzschild anti-de Sitter space referenced 
to anti-de Sitter space and de Sitter space referenced to Minkowski space. 



5 Application to black hole thermodynamics 



One application of our expressions is to black hole thermodynamics [10,12,13]. 
For this purpose we want to control the quasilocal energy-momentum so we 
must allow N to vary. Hence we Legendre transform to the "microcanonical" 
Hamiltonian iJ micro (A) := H(N) — § dTi B(N). We choose the connection as 



one of our control variables, and use (6) but drop the unphysical dynamic 
reference frame contribution due to £/v^ Q . We take our boundaries at oo and 
on the bifurcate Killing horizon. For the evolution vector field N we use the 
Killing field \ := &t + ^ud^ which is normal to and vanishes on the horizon. 
We obtain the "first law" for a general gravity theory by evaluating "on shell" 



= 5H micro (x) = j> SB( X ) = / SB( X ) - j> SB( X ) 

oo H 

(5E + n H 5J)- j>D aX ^ Pp a . (13) 



as oo H 



H 



The last integral, which for Einstein's theory has the form ~ k,5A, identifies 
the entropy for these general geometric gravity theories, 



6 Spinor expressions 



Some of the proposed quasilocal expressions for Einstein's theory are formu- 
lated in terms of auxiliary spinor fields. Here we indicate the relationship 
between such expressions and our expressions and briefly consider the appli- 
cation of our formalism to spinor formulations. 

Via certain new spinor- curvature identities [14] several new quadratic spinor 
Lagrangians for Einstein's theory have been found [15]. The different versions 
depend on whether the connection is varied independently and how the van- 
ishing torsion constraint is imposed. One of the simplest is 

£ qB :=2D(^?)7 5 lW)= -W>R*1 + 4^(^)75^+^75^ (#)}■ (14) 



It differs from the usual Einstein-Hilbert action by just a total differential. The 
variables are a Dirac matrix valued orthonormal frame one-form 1? := , ~f a 'd a 
and a "normalized" spinor field ip (i.e., ipifj = 1, ip'j^ip = 0). Asymptotically 
ip ~ const + 0(l/r) so the Lagrangian is 0(l/r 4 ) which guarantees finite 
action. 

The corresponding covariant (we have dropped a term proportional to i^oj 
which generates the frame gauge transformations and vanishes on shell) Hamil- 
tonian 3-form has the form 



H qs (N) :=2{D(W)>y 5 D(foJ>) + D(^) l5 D(^)} (15) 

= -2^/^GV * r + 2d{V# 75 £W) + £>(^)75W, (16) 
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which is just the ADM Hamiltonian up to a total differential. It is asymptoti- 
cally 0(l/r 4 ) and its variation has an 0(l/r 3 ) boundary term which vanishes 
asymptotically so there is no need for a further adjustment [3] by an addi- 
tional boundary term. Expression (15) is similar to the Hamiltonian 3-form 
associated with the Witten positive energy proof [16]: 

W w (V0 := 2(£ty 75 £W) + D(fr&) l5 Ihl>) (17) 

= -2JV"G^ *#" + 2d($fry 5 Drl> - D^ty), (18) 

wherein iV M = ifj^^ifj. Again this is the ADM Hamiltonian up to a total dif- 
ferential but TC W (unlike 7i qs ) is not related to a Lagrangian. For these spinor 
Hamiltonians once again on a solution only the boundary terms contribute to 
the value. 

Note that for these spinor expressions there is no explicit need for a reference 
configuration, the spinor field implicitly plays this role [17]. In order to com- 
pare these spinor expressions with the quasilocal expressions discussed earlier 
we introduce an explicit reference configuration. We then find that both spinor 
expressions are related to our expression (9): 

B qs (N) : = 2(# 75 5(#) + £>(^)75W) (19) 

= Au a p i N eJ - 2(^ 7 5#£)V> + I%5W), (20) 

£ W (V>) := -2(^7 5 tfDV> + £>?75#) (21) 

= Au a pi N eJ - 2(^75^ + £V75^)- (22) 

In the limit r — > 00 these spinor boundary expressions also give the correct 
total energy-momentum. Expressions like Eq. (21) have been used in several 
quasilocal energy investigations [18]. 

Having introduced a reference configuration, instead of just comparing the 
Hamiltonian boundary expression we can apply our general formalism. The 
momenta conjugate to the spinor fields can be introduced; then, in addition to 
the frame and connection type terms we had earlier, the quasilocal expression 
acquires spinor terms of the sort (3, 4) and the variation of the Hamiltonian 
contains extra spinor field contributions of the type (7). 



7 Discussion 



Which quasilocal expression gives the correct physics? The physical role of the 
spinor field especially still seems mysterious. One way to investigate these vari- 
ous quasilocal expressions is to do more direct calculations for exact solutions, 



e.g., [19]. However, a deeper theoretical investigation could be more reveal- 
ing. Our formulation provides a good starting point for such an investigation. 
Note that all of the expressions presented here correspond to the work done 
in some (ideal) physical process. The situation is similar to thermodynamics 
with its different energies (enthalpy, Gibbs, Helmholtz, etc.) An even better 
analogy is the electrostatic work required while controlling the potential on 
the boundary of a region vs. that required while controlling the charge den- 
sity [2]. Thus for the spinor expressions we simply have a different boundary 
symplectic structure corresponding to different control variables. What must 
be held fixed in each case is found by calculating the boundary term in the 
variation of the Hamiltonian. Mathematically this is straightforward. But no 
matter which technical procedure is used for the relation between the frame, 
the connection and the vanishing torsion condition the complete results for the 
variational symplectic structure turns out to be rather complicated. Briefly, 
for the spinor expressions the main conclusion is, not surprisingly, that we 
must hold the orthonormal frame and the spinor field fixed on the boundary. 
Hence to understand the physics of the quasilocal spinor expressions one must 
understand the physical meaning of controlling the spinor field on the bound- 
ary. For the Witten Hamiltonian the relation N^ = i/jj^i/j already gives part 
of the answer. 

Thus, for general geometric gravity theories, from a covariant Hamiltonian 
formulation using differential forms, by always working with four dimension- 
ally covariant variables and by using symplectic ideas, we have found several 
manifestly covariant Hamiltonian boundary term expressions for the quasilo- 
cal quantities: energy-momentum and angular momentum. Our quasilocal ex- 
pressions depend only on observer independent geometric quantities. They are 
differential 2-form expressions which can be evaluated on any closed 2-surface. 
These quasilocal expressions depend on (i) a field configuration, (ii) a refer- 
ence configuration (or spinor field), and (iii) an evolution vector field N on the 
boundary. Our formulation should be a good basis for further investigations 
aimed at understanding both the physical meaning of the quasilocal quan- 
tities for the different control modes and the physical interpretation of the 
contributions from auxiliary spinor fields. 



Acknowledgment 



We would like to thank V. V. Zhytnikov for his helpful discussions. This work 
was supported by the National Science Council of the Republic of China under 
contracts NSC 83-0208-M-008-014, 84-21 12-M-008-004. 



References 

[I] J. M. Nester, Mod. Phys. Lett. A 6, 2655 (1991); 

in: Directions in General Relativity, Vol.1, eds. B. L. Hu, M. P. Ryan and C. V. 
Vishveshwara, (Cambridge, 1993) p. 245. 

[2] J. Kijowski and W. M. Tulczyjew, A Symplectic Framework for Field Theories 
(Lecture Notes in Physics 107, Springer- Verlag, Berlin, 1979); 

J. Jezierski and J. Kijowski, Gen. Rel. Grav. 22 (1990) 1283; 

J. Lee and R. M. Wald, J. Math. Phys. 31 (1990) 725; 

D. Sudarsky and R. M. Wald, Phys. Rev. D 46 (1992) 1453. 

[3] T. Regge and C. Teitelboim, Ann. Phys. 88 (1974) 286. 

[4] R. Hecht, Mass and Spin of Poincare Gauge Theory, Gen. Rel. Grav. (1995), to 
be published. 

[5] C. M. Chen, Quasilocal Quantities for Gravity Theories, MSc. Thesis, National 
Central University, 1994. 

[6] J. Kijowski, in: Proceedings of Journees Relativistes 1983, Torino, ed. S. Benenti 
et al. (Pitagora Editrice Bologna, 1985), p. 205. 

[7] R. Beig and N. 6 Murchadha, Ann. Phys. 174 (1987) 463. 

[8] K. Hayashi and T. Shirafuji, Prog. Theor. Phys. 73 (1985) 54; 
T. Kawai, Prog. Theor. Phys. 79 (1988) 920; 
M. Blagojevic and M. Vasilic, Class. Quant. Grav. 5 (1988) 1241. 

[9] R. D. Hecht and J. M. Nester, Phys. Lett. A 180 (1993) 324; 

Evaluation of PGT mass and spin at null infinity, National Central University 
preprint, 1995. 

[10] J. D. Brown and J. W. York, Phys. Rev. D 47 (1993) 1407, 1420. 

[II] J. Katz, Class. Quant. Grav. 2 (1985) 423; 

J. Katz and A. Ori, Class. Quant. Grav. 7 (1990) 787. 

[12] J. D. Brown, J. Creighton and R. B. Mann, Phys. Rev. D 50 (1994) 6394. 

[13] R. M. Wald, Phys. Rev. D 48 (1993) 3427; 

T. Jacobson and R. C. Meyers, Phys. Rev. Lett. 70 (1993) 3684; 

T. Jacobson, G. Kang and R. C. Myers, Phys. Rev. D 49 (1994) 6587. 

[14] J. M. Nester, R. S. Tung and V. V. Zhytnikov, Class. Quant. Grav. 11 (1994) 
983. 

[15] J. M. Nester and R. S. Tung, Gen. Rel. Grav. 27 (1995) 115. 



[16] J. M. Nester, in: Asymptotic Behavior of Mass and Space-time Geometry, ed 

F. Flaherty (Springer, Berlin, 1984) p. 155. 

[17] S. Lau, Class. Quant. Grav. 12 (1995) 115. 

[18] A. J. Dougan and L. J. Mason, Phys. Rev. Lett. 67 (1991) 2119; 

G. Bergqvist, Class. Quant. Grav. 9 (1992) 1753, 1917. 

[19] E. A. Martinez, Phys. Rev. D 50 (1994) 4920. 



10 



